Properties of J/tp at T c : QCD second-order Stark effect 
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Starting from the temperature dependencies of the energy density and pressure from lattice QCD 
calculation, we extract the temperature dependencies of the electric and magnetic condensate near 
T c . While the magnetic condensate hardly changes across T c , we find that the electric condensate 
increases abruptly above T c . This induces a small but an equally abrupt decrease in the mass of 
J/tp, which can be calculated through the second-order Stark effect. Combining the present result 
with the previously determined QCD sum rule constraint, we extract the thermal width of J/tp 
above T c , which also increases fast. These changes can be identified as the critical behavior of J/ if) 
across T c associated with the phase transition. We find that the mass shift and width broadening 
of J/tp at 1.05 T c will be around -100 MeV and 100 MeV respectively. 
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Since the pioneering work by Hashimoto et al. [l[ and 
the seminal work by Matsui and Satz (2j, many ex- 
perimental and theoretical works have been performed 
in the physics of J/tp suppression in heavy ion colli- 
sions. The subject has recently evolved into a new stage 
as lattice calculations using maximum entropy methods 
(MEM) found the peak structure to survive up to al- 
most 2T C [H, H| , which was speculated before as lat- 
tice calculations showed that non-perturbative nature of 
QCD persists well above T c [f| 0j. This suggests that 
the sudden disappearance of J/tp is not the direct sig- 
nature of QGP formation. Indeed, recent results from 
RHIC on the suppression factors at different rapidity and 
higher px seem ever more confusing [§[ . Hydrodynamic 
calculations suggest that the initial temperature of the 
QGP formed at RHIC is in the order of 2T C and will last 
for 3 to 4 fm/c @, [l(J. Therefore, considering the for- 
mation time of charmonium states after the creation of 
cc pair, it is crucial to know the detailed properties of 
J/tp near T c to fully understand the suppression and/or 
enhancement of J/tp in heavy ion collisions. Unfortu- 
nately, the present resolutions of the peak structure of 
J/tp from the lattice calculations based on MEM are far 
from satisfactory In fact, the peak is too broad 

to even discriminate between J/tp from tp' . Moreover, 
the temperature region between T c to 2T C is known to 
be strongly interacting and therefore a non-perturbative 
method has to be implemented to consistently treat the 
charmonium at this temperature region. In a previous 
work [T3]j we have implemented QCD sum rules to in- 
vestigate the properties of J/tp near T c . Although the 
results were non-perturbative, only a constraint on the 
combined mass decrease and width increase could be ob- 
tained. Here, we point out that the critical behavior of 
QCD phase transition, could be identified with a critical 
behavior of electric condensate at T c , and then by making 
use of the QCD second order Stark effect, show that such 
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critical behavior can be translated to a sudden change in 
the mass of J/tp across T c . 

We begin by characterizing the properties of the 
strongly interacting quark-gluon plasma (sQGP) of the 
pure gluon theory across T c in terms of local opera- 
tors. This is accomplished by making use of the energy- 
momentum tensor, which has a symmetric traceless part 
and a trace part via the trace anomaly, 



T. 



a/3 



-ST(G^GJ) 



9a/3 P(g) 



Here, a and a, p3 are the color and Lorentz indices re- 
spectively. The temperature dependence of the two in- 
dependent parts of the energy-momentum tensor can be 
obtained from the lattice measurement of energy density 
and pressure at finite temperature. 



(T af3 ) T = (e+p) 



1 



9al3 + (e - 3p) 



9af3 



(2) 



Here, u a is the four velocity of the heat bath. Therefore 
the temperature dependencies of gluonic operators can 
be identified with the pressure and energy density. To 
leading order in coupling, we can identify the trace part 
-f [<V G ^ Ga HT - (fG°„G a no] = M (T) and the 
non-trace part {~STG a a ^G a ^)) T = (u a u p -\g a p)M 2 (T), 
where [13, 



M (T) = (e-3p), 
M 2 (T) = (e+p), 



(3) 



and (v G ^ GaMiy )o = G o ac is thc scalar S luon conden- 
sate in the vacuum. The lattice gauge theory result 
for e and p in the pure SU(3) gauge theory was ob- 
tained from Rcf. hj. Fi gurc Q] shows changes of Mq 
and Mi from their vacuum value scaled by their asymp- 
totic temperature dependence of T 4 . One notes that 
while M-2, which is also proportional to entropy density 
times temperature, moderately reaches the asymptotic 
temperature dependence, Mq, also known as the interac- 
tion measure or the gluon condensate, suddenly increases 
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FIG. 1: M = (e - 3p) and M 2 = (e + p) divided by T 4 as 
functions of temperature. 



FIG. 2: Electric and magnetic condensate near T c as functions 
of temperature. 



and then decreases at higher temperature. The strongly 
interacting nature of QGP is related to the large inter- 
action measure (l5j . which takes its maximum value at 
around 1.1T C . It should be noted that the temperature 
dependence of the gluon condensate M extracted here 
includes both the perturbative and the non-perturbativc 
and thus the full temperature dependence. The sudden 
change near T c is dominated by the sudden decrease of 
the non-perturbativepart, which reduces to about half 
of its vacuum value [7(, while at higher temperatures it 
is dominated by the perturbative contributions [H, [lj| . 

For the heat bath at rest, one can rewrite the ther- 
mal expectation values of the dimension four operators 
of the energy-momentum tensor in Eq. ((T|) in terms of 
electric and magnetic condensate [l7[. This is possible 
after making the following identification. 



(4) 



The scale dependence of the matrix element is transferred 
to the coupling constant. Therefore, we additionally need 
to know the temperature dependence of the coupling con- 
stant a s (T). Since we will be using the matrix element in 
the operator product expansion (OPE) with the separa- 
tion scale relevant for the heavy bound state, we will 
use the temperature dependent running coupling con- 
stant extracted from the lattice computation of the heavy 
quark free energy Then we find, 
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Figure [5] shows the temperature dependence of 
{^E 2 ) T and (^B 2 ) T - One notes that there is a sudden 



increase in the electric condensate ( 9 ^-E 2 )t, while the 
magnetic condensate {^-B 2 )t hardly changes above T c . 
This can be related to the fact that the area law behavior 
of the space-time Wilson loop changes to the perimeter 
law above T c , while that of the space-space Wilson loop 
retains the area law behavior even above T c Q . The con- 
nection comes in as the non perturbative behavior of a 
rectangular Wilson loop in the SiS% direction can be re- 
lated to the non- vanishing gluon condensate (^G| iS2 ) 
via the operator product expansion Hence, one can 
conclude that critical behaviors of QCD phase transi- 
tion can be related to the sudden change in the electric 
condensate. Such local changes will induce critical be- 
havior of a heavy quark system such as the J/ip across 
the phase transition, which can be obtained through the 
QCD second-order Stark effect. 

The perturbative QCD formalism for calculating the 
interaction between heavy quarkonium and partons was 
first developed by Peskin [lj|, [2(| in the non-relativistic 
limit. The formula for the mass shift reduces to the 
second-order Stark effect in QCD, which was used pre- 
viously to calculate the mass shift of charmonium in 
nuclear matter [l7j . The information needed from the 
medium is the electric field square. As the dominant 
change across the phase transition is the electric conden- 
sate, one notes that the second-order Stark effect is the 
most natural formula to be used across the phase transi- 
tion. 

The second-order Stark effect for the ground state 
charmonium with momentum space wave function nor- 

p) I 2 = 1 is as follows, 
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FIG. 3: Mass shift from the second-order Stark effect (solid 
line) and the maximal mass shift obtained from QCD sum 
rules from Ref. [l2| (points). 



where k = \k\ and {^-AE 2 )t denotes the value of change 
of the electric condensate from its vacuum value. The sec- 
ond line is obtained for the Coulomb wave function. Here, 
e is the binding energy and m c the charm quark mass. 
These parameters are fit to the size of the wave function 
obtained in the Cornell potential model [2l|, and to the 
mass of J/tjj assuming it to be a Coulombic bound state 
in the heavy quark limit [l9[. The fit gives m c = 1704 
MeV, a = 0.271 fm and a s = 0.57. Few comments are in 
order. The minus sign in Eq. |(7J) is a model independent 
result and follows from the fact that the second-order 
Stark effect is negative for the ground state. The factor 
of Bohr radius square a 2 follows from the dipole nature of 
the interaction, and the binding energy e from the inverse 
propagator, characterizing the separation scale [l9l. I22I]. 
Therefore, the actual value of the mass shift does not de- 
pend much on the form of the wave function as long as 
the size of the wave function is fixed. The Bohr radius 
used in our calculation corresponds to (r 2 ) 1 / 2 = 0.47 fm, 
which is the size of a more realistic wave function in the 
Cornell potential [2l[. Therefore, the correction coming 
from using a more realistic wave function should be small. 

The solid line in Fig. [3J shows the mass shift obtained 
from the second-order Stark effect. The second-order 
Stark effect formula is based on the operator product 
expansion (OPE) for bound state. As mentioned before, 
the formalism was first established by Pcskin in 1979 and 
the separation scale is the binding e = mg A . A more sys- 
tematic derivation was developed recently by Brambila 
ct al. [23| for the bound state. Here, the relevant scales 
are mv and mv 2 , where the former is related to the po- 
tential 1/r and later the kinetic energy p 2 jm. This scale 
mv 2 is the separation scale so that for effects with typ- 
ical momentum larger than the separation scale should 
be taken into account through resummed perturbation 
by solving the Schrodinger equation (22j, while that with 
smaller scale should be taken into account through the 



operator product expansion. The operators contain the 
non-perturbative physics of QCD, which is typically of or- 
der Aqcd ■ Therefore the OPE for the bound state works 
best when mv 2 3> Aqcd [HI- While there are concerns 
that this condition is marginal for charmonium, the ap- 
proach should provide a quantitative description. Now 
the question is, which approach should one take for the 
thermal interactions near T c . Obviously, the effects of fi- 
nite temperature involves new scales like the color screen- 
ing. However, as has been known for some time, the tem- 
perature region from T c to 2 T c , is known to be strongly 
interaction and can not be described by resumed pertur- 
bation [24j| . Therefore, the approach we want to take is 
to calculate the non perturbativc temperature effect to 
the mass shift through the operator product expansion. 
The leading order contribution in this approach is coming 
from (^AE 2 ) T as is given in Eq. ([7]), whose temper- 
ature dependence we extract directly from the lattice. 
The arguments for convergence of higher dimensional 
operators in our approach are twofold. First, we will 
restrict to the temperature region where the change in 
(^AE 2 ) T is smaller than the vacuum value of (^E 2 ) Q 
itself. We believe that then the OPE is under control 
as has been verified by the typical QCD sum rule ap- 
proaches for heavy quark system in the vacuum. As can 
be seen in Fig. [2j this condition restricts our applicabil- 
ity to 1.05T C , as in the QCD sum rule approach at finite 
temperature [12, HH|. Second, a more direct evidence 
comes from the next term in the OPE correction, which 
comes from magnetic condensate. However, as can be 
seen from Fig. [21 the changes of ^AB 2 ) T and hence 
the next term in the OPE should be small up to 1.05T C 
and slightly beyond. Therefore, we can conclude that the 
second-order Stark effect should be valid near T c . As can 
be seen in Fig. [31 the results from second-order Stark ef- 
fect shows that the mass reduces abruptly above T c and 
becomes smaller by about 100 MeV at 1.05T C , reflecting 
the critical behavior of the QCD phase transition. 

We put the present result in perspective with a non- 
perturbative result obtained before using the QCD sum 
rules (l2l . I25I ]. The points in Fig. [3j represent the max- 
imum mass shift obtained in Refs. [l2|, [25J. As can 
be seen in the figure, the mass shift obtained from the 
second-order Stark effect is almost the same as the max- 
imum mass shift obtained in the sum rule up to T c 
and then becomes smaller. The mass shift at T c is 
about —50 MeV. In the QCD sum rules, only a con- 
straint for the combined mass shift and thermal width 
of -Am + T T ~ 80 + 17(T - T c ) MeV could be ob- 
tained within the temperature range from T c to 1.05T C . 
Therefore, the difference between the Stark effect and the 
maximum mass shift obtained from QCD sum rules above 
T c in Fig. [3j could be attributed to the non-perturbative 
thermal width at finite temperature. In Fig. [H we plot 
the thermal width obtained from combining the QCD 
sum rule constraint with the mass shift obtained from 
the QCD second-order Stark effect. [37| As can be seen in 
the figure, the thermal width at 1.05T C becomes as large 
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FIG. 4: Thermal width of J/^/j obtained from the second-order 
Stark effect and QCD sum rule constraint. 



as 100 MeV. Such width slightly above T c is larger than 
that estimated from a perturbative LO and NLO QCD 
method [26|, [27j , but smaller than a recent phenomeno- 
logical estimate J28|. The mass of quarkonium at finite 
temperature was also investigated in the potential mod- 
els [29|], where the mass was found to decrease at high 
temperature. However, the detailed potential has to be 
extracted from the lattice 1 at each temperature and hence 
identifying the critical behavior near T c will be difficult. 

Finally, we comment that while the present result is 
obtained using lattice calculation in the pure gauge the- 
ory, including dynamical quark will not greatly modify 
the result. This follows from noting that the lattice re- 
sult for the temperature dependence of pressure p/T 4 is 
independent of the number of flavors if scaled to their cor- 
responding ideal gas limit (30l |. Moreover, as was shown 
in Ref. [31[ , Go (T) extracted from a recent full lattice cal- 
culation of the interaction measure [32| after subtracting 
the quark contributions, and then dividing by a factor of 



(l + j^rif) appearing in the beta function, shows that the 
change of the magnitude near T c is remarkably similar to 
that of the pure gauge theory. Hence the main input for 
our result does not change much even in the presence of 
dynamical quarks. 

In summary, we have shown that the sudden increase 
in the energy density across the phase transition, which 
is a characteristic behavior of the QCD phase transition 
independent of the flavor, can be translated to a rapid in- 
crease in the electric condensate slightly above T c . Using 
the QCD second-order Stark effect, this translates into 
an equally sudden decrease in the mass of J/ip, which 
is around —50 MeV and —100 MeV respectively at T c 
and 1.05T C . Combining with a QCD sum rule constraint, 
we obtain the thermal width of J/ip slightly above T c , 
and found it to be larger than previous perturbative es- 
timates, and becoming as much as 100 MeV at 1.05T C . 
Hence, one can conclude that the critical behavior of J/ip 
at T c is not its sudden disappearance, but rather the 
abrupt changes of its mass and width. The mass shift 
is probably too small to be detected with present reso- 
lutions at RHIC. However, with the expected upgrades 
at RHIC and plans at LHC, such direct measurement 
could be possible. Indeed, the mass resolution of J/ip 
for dimuon channel at LHC is 35 MeV for the CMS de- 
tectorjp and around 70 MeV for ALICE [H[ and AT- 
LAS [35[ . It might be better for diclectron channel. Fur- 
thermore, the mass shift could also influence production 
rates within the statistical model [36| . The large width 
of 100-150 MeV already at 1.05T C suggest that while the 
maximal entropy method shows a J/ip peak structure 
surviving up to 2T C , the actual formation at heavy ion 
collisions might only be possible at lower temperatures 
where the width of the J/ip becomes equal to its bind- 
ing. 
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